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FOREWORD 


The  analyses  described  in  this  report  are  based  on  a  problem  originated 
by  Mr.  Said  Saadi  (W60),  namely  to  develop  procedures  for  filling  in  data 
gaps  in  the  MARS  (Mission  Analysis  Reporting  Suite)  data  base.  The  study 
carried  out  is  exploratory  in  nature  but  does  give  a  beginning  to  addressing 
Mr.  Saadi’s  problem. 

The  author  is  indebted  to  Mr.  Saadi  for  helpful  discussions  and  a  review 
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I.  INTRODUCTION 

The  technology  for  intercepting  a  threat  missile  (TM)  by  an  interceptor 
missile  (IM)  requires  numerous  ground  and  flight  tests  to  establish  its  valid¬ 
ity.  A  huge  amount  of  information  is  compiled  from  these  tests  and  is,  in  a 
particular  case,  referred  to  as  the  MARS  (Mission  Analysis  Reporting  Suite) 
data  base.  Often  ’’holes”  are  observed  in  the  data,  that  is,  some  desired  data 
is  not  contained  in  the  data  base.  The  objective  in  this  report  is  to  give  an 
example  of  how  an  interpolation  procedure  can  supply  estimates  for  missing 
data  using  only  the  available  experimental  data. 

The  example  chosen  is  based  on  simulated  notational  data  rather  than  on 
the  MARS  data,  since  the  latter  is  classified.  However,  the  simulation  is 
designed  to  be  more  of  a  challenge  for  interpolations,  since  the  data  from  the 
simulation  are  less  smooth  than  that  generated  from  MARS. 

The  example  is  taken  from  the  following  scenario,  with  all  action  taking 
place  exoatmospherically,  unless  noted  otherwise:  A  threat  missile  (TM)  is 
launched  from  Asia  with  Washington,  D.C.  as  the  target.  Once  it  is  acquired, 
time  is  set  to  zero  and  an  IM  is  launched  from  Alaska.  It  is  assumed  to  fly 
in  powered  flight  from  launch  for  tpow  seconds,  say  tpow  =  158  seconds. 
At  that  time  the  IM  will  have  arrived  at  a  burnout  position  (BP)  with  a 
given  velocity  (BV).  These  quantities  are  determined  by  a  three-dimensional 
flyout  fan  of  trajectories,  usually  supplied  by  the  IM  designer.  A  particular 
trajectory  is  chosen  from  the  flyout  fan  giving  a  BP  that  depends  on  the  time 
of  flight  for  the  IM  to  fly  ballistically  to  the  threat  at  position  (TP).  Then 
BV  is  determined  by  an  iterative  process  to  result  in  intercept.  Thus,  at  a 
time  ti  beyond  tpow,  when  intercept  is  first  possible,  a  “time  launch  window” 
(TLW)  is  opened  and  continues  to  stay  open  as  time  increases,  until  intercept 
is  no  longer  feasible,  say  at  tc  when  the  threat  reenters  the  atmosphere  or 
perhaps  the  BP  and/or  BV  required  at  that  time  can  no  longer  be  achieved. 
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It  is  assumed  then  that  it  can  be  first  intercepted  at  the  time  of  flight,  t  =  1 1, 
and  continuously  at  all  times  thereafter  until  it  is  no  longer  possible  at  tc. 
The  TLW  is  defined  by  the  time  interval  [ti,tc]  . 

The  problem  specifically  set  then  is  to  develop  vector  time  functions,  de¬ 
pendent  on  the  time  to  intercept,  threat  position,  threat  velocity,  and  closing 
velocity  of  the  IM,  that  can  be  used  as  interpolation  functions  for  the  burnout 
positions,  BP,  and  burnout  velocities,  BV,  of  the  IM  at  any  specified  times 
in  the  TLW. 

II.  INTERPOLATION  FUNCTIONS 

The  interpolation  functions  will  require  estimating  some  linear  and  some 
nonlinear  parameters.  The  linear  parameters  will  be  determined  by  least 
squares  (LS)  and  the  nonlinear  ones  by  assigning  them  values  over  a  spectrum 
of  possibilities  to  be  described. 

At  M  times,  in  equal  increments,  At,  throughout  the  TLW,  M  BP  ’s  and 
M  BV ’s  required  for  intercept  are  assumed  available  as  input  for  [TO]  = 
(ti,  ti  +  At,  ti  +  2At . . . ,  ti  +  (M  —  2)At,  tc).  We  will  also  use  the  notation 

[TO]  =  [t (k)] ,  t(k)  =  ti  +  (k  —  l)At,  t(M)  =  tc  k  =  1,  2, . . . ,  M  —  1,  (1) 
where  the  [.]  brackets  are  used  to  denote  a  column  vector.  The  equal  time 
increment  in  [TO]  is  used  to  deal  with  a  desirable  uniform  collection  of  data, 
but  it  is  not  a  critical  requirement. 

Additional  inputs  required  are  the  TP,  the  threat  velocity,  TV,  and  the 
closing  velocity  of  the  IM,  CV,  for  each  t  in  TO.  Thus,  for  definiteness,  we 
seek  estimates  for  the  burnout  position,  BP  =  (bx,  by,  bz),  and  velocity  of 
the  IM,  BV  =  (bu,  bv,  bw),  for  t  in  [ti,  tc]  which  includes  [TO].  The  argument 
(k)  appears;  it  always  refers  to  the  kth  time  of  [TO].  For  example,  t(k)  refers 
to  the  kth  element  of  [TO] . 

Additional  notation  is  introduced  to  simplify  expressing  subsequent  math¬ 
ematical  quantities.  Let 
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[bX(k)] 

= 

bx 

[b2  (k)] 

— 

by 

[b3  (k)] 

— 

bz 

[b4(k)] 

— 

bu 

[b5  (k)] 

— 

bv 

[b6  (k)] 

— 

bw 

[xi  (k)] 

— 

[TO] 

[x2  (k)] 

— 

[tx] 

[x3  (k)] 

— 

[ky] 

[x4(k)] 

— 

[tz] 

[x5  (k)] 

— 

[tu] 

[x6  (k)] 

— 

[tv] 

[M(k)] 

— 

[tw] 

[x8  (k)] 

— 

H 

[x9  (k)] 

— 

[cv] 

[xio(k)] 

— 

[cw] 

IM  kth  x  —  position  coordinate  at  burnout 
IM  kth  y  —  position  coordinate  at  burnout 
IM  kth  z  —  position  coordinate  at  burnout 
IM  kth  x  —  velocity  component  at  burnout 
IM  kth  y  —  velocity  component  at  burnout 
IM  kth  z  —  velocity  component  at  burnout 
M,  IM  times  of  flight  in  TLW  to  intercept 
M,  x  —  position  coordinates  of  the  threat 
M,  y  —  position  coordinates  of  the  threat 
M,  z  —  position  coordinates  of  the  threat 
M,  x  —  velocity  components  of  the  threat 
M,  y  —  velocity  components  of  the  threat 
M,  z  —  velocity  components  of  the  threat 
M,  closing  x  —  velocity  components  of  the  IM 
M,  closing  y  —  velocity  components  of  the  IM 
M,  closing  z  —  velocity  components  of  the  IM. 


(2) 


Each  [xj(k)]  is  a  column  array  with  M  specified  elements  and 
Xj  =  mean([xj(k)]).  For  example:  At=100,  [xi(k)]  =  [900, 1000, 1100, 1200, 
...,2520],  M  =  18,  xi  =  1745.556.  Note:  TLW  =  [900,  2520], 

The  Euclidean  norm  of  the  vector  j3  with  M  real  components  is  defined  by 


£/5002 


(3) 


Lk=l 


Now  instead  of  using  the  Xj  in  the  LS  equations  below,  we  have,  using  the 
mean  and  the  Euclidean  norm  for  each  Xj,  the  new  normalized  dimensionless 
variables: 


Xi(k)  =  [xi(k)  —  xj/||  [xi  —  xj || ,  i  =  1, ...,  N,  N  =  10, 
Bm(k)  =  [bm  (k)  —  bm]  / 1 1  [bm  —  bm]  1 1 ,  m  =  1, ... ,  6,  k=l,...,M. 
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The  argument  (k)  of  a  function  refers  to  the  kth  element  of  [TO] .  It  is  assumed 
that  M  data  points  of  all  the  normalized  variables  of  (4)  are  available.  The 
LS  analysis  allows  nonlinear  expressions  in  these  variables.  Nevertheless,  the 
problem  of  determining  the  LS  constants  remains  linear  and  is  classical. 

The  LS  equations  follow,  with  i  =  1, . . . ,  N  and  j  =  i, . . . ,  N.  Also  two 
integer  variables  Jm,  Jlm  are  introduced,  with  Jm  e  [1,  N —  1] ,  Jlm  e  [Jm+1,  N)] 
together  with  a  variable  T“  that  can  take  three  possible  values,  for  each  m. 
They  are  used  as  exponential  fitting  parameters  with 


7im 

if 

1  ^  n  <  Jm 

if 

if 

Jm  <  n  ^  Jlm 

(5) 

if 

if 

Jlm  <  n  ^  N,  m  =  1, ...,  6. 

The  role  of  Jm  and  Jlm  will  be  described  below.  The  interpolation  time 
functions  for  Bm  are  determined  by  the  RH  side  of  equations  (6)- (12). 


Si,  j  (k) 
B1  (k) 

B2(k) 

B3(k) 

B4(k) 

B5(k) 

B6(k) 


=  sign(Xi(k)*Xj(k)),  k  —  1, . . . ,  M 

N  N  N 

~  E * x*(k>  +  E E Ai.i  Suck)  ixi(k)ir.1  ixj(k)|U 

i=l  i=l  j=i 

N  N  N 

~  Eb‘X>(k)  •  VVlViS..:ki  |Xs(k)|rf  |Xj(k)|r? 

i=l  i=l  j=i 

N  N  N 

~  E^W  +  EEc'isy<k>  tx>(k)lr?  txi(k)lrf 

i=l  i=l  j=i 

N  N  N 

~  Ed‘x‘(k)  +  EE°y si.j(k) lxi(k)lr? lxj(k)|rf 

i=l  i=l  j=i 

N  N  N 

~  Ee>X‘(k)  +  E  EE«  Si, j (k)  |Xi(k)jri  |Xj(k)|rf 

i=l  i=l  j=i 

N  N  N 

~  Ef>x><k)  +  EEFus>.i(k)  ixi(k)ir?iXj(k)ir'. 

i=l  i=l  .1=  i 


(6) 

(7) 

(8) 
(9) 

(10) 

(11) 

(12) 


4 


NSWCDD/TR— 10/113 


The  system  of  equations  given  by  (7)- (12)  can  be  written  as  one  expression, 
i.e., 

N  N  N 

Bm(k)  ~  £  aT  X,(k)  +  J2  A”  Z”  00.  (13) 

i=l  i=l  j=i 

where 


ai  —  ai)  ai  —  aj  —  fi,  Ay  —  Ay,  Ay  —  By,  ...,  Ay  —  Fy , 

Z“j(k)=SiJ(k)|Xi(k)|r"  |Xj(k)|pj”. 


(14) 

(15) 


Thus,  assigning  a  value  in  (13)  to  m,  say  3,  would  also  refer  explicitly  to  (9). 
Also,  to  clarify  the  use  of  Jm  and  Jlm  with  T“,  (5),  say  for  m=l,  the  following 
quantities  appear  in  (7): 


IX/S  |X2f’,  iXal1-,  |X4|i'S  \Xs\L\  | X6 1 1  ® ,  |X7|^,  iXg^S,  l^l^,  IX^X  (16) 


iri 


iri 


iri 


\Tl 


r1 

-L  R 


irf 


iri 


iri 


irio 


Then,  for  example,  if  J1  =  3,  and  Jl1  =  7,  (16)  becomes,  from  (5), 

iXrPs  |X2rs  |X3rs  |X4p\  |X5p\  |X6p\  |X7p\  iXgp1,  |X9p\  |X10p».  (17) 


From  (13)  we  see,  for  each  m,  there  are  NT  =  N  +  N(N  +  l)/2  linear 
parameters  to  determine  and,  in  addition,  five  others,  namely, 


jrm 


(Jm,  Jlm,7im,72m,73m), 


(18) 


for  a  total  of  70  when  N=10.  Assuming  the  elements  of  JT111  are  assigned 
values,  then  for  each  value  of  m  =  1, ...,  6  there  is  an  LS  problem  to  solve  for 
the  variables 


dl  ,  d2  j 


(NTx  1) 


111  A  111  A  111 

aN  i  "‘M,!  ’  ^1,2) 


A  111  A  111  A  111 

a1,N>  a2,2  5  ^2,3? 


(19) 


The  right  hand  side  of  the  LS  equations  is  given  by 

(Mxl) 

Bm  =  [Bm(l),  Bm(2), . . . ,  Bm(M)],  (20) 
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and  the  LS  matrix  is  given  by 

(MxNT) 

|"  Xi  (1)  ...  XN(1)  Zft(l)  ...  Z^(l)  ...  Z^(l)' 

Xi(2)  ...  XN(2)  Zft(2)  ...  Z™n(2)  ...  Z^(2) 

Xi(M)  ...  XN(M)  Zft(M)  ...  Z^N(M)  ...  ZJ>N(M) 

Hence,  in  matrix-vector  form,  we  have 

XZm£m  =  Bm  (22) 

Solving  this  system  of  linear  equations,  one  obtains 

C  =  (XZ“  XZm)~1XZ“  Bm,  (23) 

where  XZ^  denotes  the  transpose  of  XZm. 

The  maximum  error  generated  by  an  LS  analysis  generally  occurs  in  the 
vicinity  of  midway  between  the  data  points.  In  this  light,  we  employ  addi¬ 
tional  input  similar  to  the  original  normalized  input,  as  given  by  (4).  This 
new  input  data  will  be  taken  at  points 

tl(l)  =  (t(l)  +  t(2))/2,  [Tl]  =  (tl(l),  ...,tl(M  —  1)), 
tl(kl)  =  tl(l)  +  (kl  -  l)At,  kl  =  1,  ...,M  -  1. 

and  is  denoted  by 

Xli(kl)  =  [xli(kl)  -Xi]/||xli  -  Xi||,  i  =  1, ...,  N  =  10, 

Blm(kl)  =  [blm(kl)  —  b™"]/||blm  —  b“j|,  m  =  1, ...,  6. 

The  argument  (kl)  of  a  function  refers  to  the  klth element  of  [Tl]. 

It  remains  to  explain  how  the  elements  of  JTm  are  chosen.  For  a  fixed  m, 
say  m  equals  one,  the  integer  variables  J1  and  Jl1  are  sequenced  through  the 
integers  1  to  N-l  and  J1  +  1  to  N,  respectively.  Then  71,72,73  are  varied 
each  in  turn  through  a  range  of  real  values,  say  [0,5],  at  initial  increments 
of  A  =  .2.  At  each  assignment  of  these  variables,  J1,  Jl1,  7J,  72,  73,  an  LS 
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solution  over  [TO]  is  carried  out  and  the  absolute  maximum  error  (AME)  over 
[TO],  ErrTO,  and  [Tl],  ErrTl,  is  noted.  This  procedure  is  iterated  varying 
J\  Jl\  7i,  72,  73,  systematically  as  described  above.  The  best  values  for  the 
elements  of  JT1  are  those  that  minimize  the  AMEs  over  the  iterations.  The 
entire  process  is  repeated  with  reduced  ranges  for  each  variable  with  a  smaller 
A.  For  example,  say,  J1  =  2,  Jl1  =  8,7!  =  2.0,  72  =  .75,  73  =  1.0;  the  ranges 
of  these  elements  are  then  reduced  by  setting  J1  =  [1,3],  Jl1  =  [7,  9] ,  7^  = 
[1.8,  2. 2],  72  =  [.55,  .95],  73  =  [0.8, 1.2]  and  A  =  A/2  =  0.1,  obtaining,  in 
general,  a  smaller  overall  AME.  The  iteration  process  is  again  carried  out 
with  A  =  .05  and  .025  with  the  final  results  after  minimizing  the  AME  given 
by,  say,  JT1  =  (J1  =  3,  Jl1  =  7, 7/  =  1.925, 7]  =  .700, 731  =  2.120). 

The  AME  over  [TO]  is  defined  as 

N  N  N 

ErrTO  =  max  |Bm(k)  —  E  a!”Xi(k)  —  EE  Z™(k)|,  (26) 

i=l  i=l  j=i 


and  the  AME  over  [Tl],  using  (24)  and  (25),  is 

N  N  N 

ErrTl  _  rnax  |B1  (kl)  E  afXli(kl)-  EEA"Zl"(kl)l, 

where 


i=i 


i=l  j=i 
rp  1  vi  ur? 


ZlPj(kl)  =  Si,j(kl)|Xli(kl)|1'  jXljfkl)!1)  . 


(27) 

(28) 


III.  NUMERICAL  EXAMPLE  WITH  GRAPHS 


Considering  the  example  described  in  the  Introduction,  we  are  now  in  a 
position  to  report  some  numerical  results.  We  start,  as  described,  with  a 
TM  flying  from  Asia  with  Washington,  D.C.  as  a  target,  with  the  IM  launch 
site  in  Alaska.  It  is  assumed  the  powered  flight  of  the  IM  takes  158  seconds 
to  reach  a  burnout  position.  A  fictitious  flyout  fan  is  used;  thus  all  data 
is  unclassified.  The  threat  position  at  t=0  is  (-1456954,5437427,3250043) 
and  its  velocity  is  (296,-1106,7404).  The  Alaska  launch  site  coordinates  are 
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(-2334397,-1347764,5780581).  All  data  refers  to  an  ECEF  frame.  Position 
coordinates  are  in  meters  and  velocity  coordinates  are  in  meters  per  second. 
The  input  data  is  generated  from  Fortran  programs  described  in  [1]  under 
Case  #1. 

The  threat  is  first  acquired  at  time  zero.  The  intercept  window  opens 
900  seconds  later  and  remains  open  until  the  window  closes  when  the  threat 
reenters  the  atmosphere  at  t=2520  seconds.  So  data  for  the  Xj  (see  (2))  will 
be  taken  in  equal  time  intervals  (see  (2),  (4)),  At  =  20  from  t=900  to  t=2520. 
Thus  M  =  82.  This  array  of  times  is  denoted  by  [TO]. 

In  the  previous  section,  two  sets  of  time  data  were  introduced,  belonging 
to  TLW,  namely,  [TO],  as  above,  and  [Tl] .  The  two  time  sets  were  made  up 
of  times  at  equal  intervals,  namely, 


CM 

i— i 

II 

o’ 

tL 

,  t]yt  —  tc, 

(29) 

[Ti]  =  (tl+t2.t2+t3,. 

1  +  I'M  \ 

2 

(30) 

Assigning  specific  values,  we  have 

[TO]  =  (900,  920,  ..., 

2500,  2520), 

(31) 

[Tl]  =  (910,  930,  ..., 

2490,  2510). 

(32) 

Associated  with  each  set  were  other  input  quantities  (normalized),  Xj(k) 
and  Xli(kl),  respectively,  for  each  time,  i.e.,  for  each  k  and  kl  (see  (4), 
(24),  (25)).  At  each  of  these  intercept  times  belonging  to  [TO]  and  [Tl],  it  is 
assumed  the  following  ECEF  data  are  available: 

IM  Burnout  position,  BP 
IM  Burnout  velocity,  BV 
Threat  position,  TP 
Threat  velocity,  TV 
IM  Closing  velocity,  CV  . 

Then,  using  this  data,  a  set  of  functions  is  developed  for  which  least  square 
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parameters  are  determined  together  with  the  assigned  nonlinear  parameters, 
jrm,  (18).  These  functions  are  then  tested  by  using  them  to  interpolate  for 
BP  and  BV  at  a  set  of  times 

3t  i  +  t2  ti  +  3t2  3t2  +  t3  3tM-l  +  tM  tM— 1  +  3tM 

4  ’  4  ’  4  ’ " '  ’  4  ’  4 

[T2]  =  (905,  915,  925,  ...,  2505,  2515),  (34) 

as  well  as  at  times  [TO]  and  [Tl] .  Corresponding  to  (25) 

X2j(k2)  =  [x2i(k2)-xi]/||x2i-xi||,  k2  =  l,...,2M-2, 

B2“(k2)  =  [b2m(k2)  -b“]/||b2m  -b“||,  m  =  1, ...,  6. 

where  x2j  is  the  same  function  as  Xj  but  evaluated  at  times  [T2]  rather  than 
[TO].  The  subscript  s  of  B2“(k2)  indicates  that  B2  comes  from  simulated 
data  in  reference  [1]  and  is  used  here  as  “truth  data.”  The  argument  (k2)  of  a 
function  refers  to  the  k2th  element  of  [T2].  The  errors  in  estimating  B2“(k2) 
are  determined  from: 

N  N  N 

ErrorT2  =  B2“(k2)  -  ^a“X2i(k2)  -  EEASZ2”(k2)'  (36) 

i=l  i=l  j=i 

Z2“(k2)  =  Sij(k2)  |X2i(k2)|r“  |X2j(k2)|rr.  (37) 

The  absolute  maximum  error  is  given  by 

N  N  N 

ErrT2  =  max  |B2“(k2)  -  a“X2i(k2)  -  EEA"iZ2"i(k2)l-  (38) 

i=l  i=l  j=i 

Programs  in  Fortran  95  and  MatLab  were  developed  to  obtain  the  results 
listed  below  in  Table  1. 
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Table  1.  Final  Values  for  JT,  (18),  and  Errors  (26),  (27),  (38) 

m 

Jm 

Jlm 

7  r 

m 

72 

m 

73 

ErrTO 

ErrTl 

ErrT2 

1 

2 

9 

4.850 

0.900 

2.275 

581 

882 

2533 

2 

5 

7 

2.275 

2.45 

1.000 

436 

784 

1348 

3 

4 

8 

1.625 

0.900 

1.000 

200 

593 

830 

4 

2 

7 

3.000 

1.500 

0.725 

0.221 

0.442 

0.820 

5 

2 

6 

3.825 

1.4000 

1.000 

0.155 

0.434 

0.590 

6 

4 

5 

2.700 

1.625 

1.000 

0.152 

0.548 

0.632 

In  Table  1,  m=2  refers  to  the  y-coordinate  burnout  position.  The  next  five 
columns  are  used  to  record  the  values  of  the  nonlinear  parameters  in  jrm  so 
that,  as  in  (16)  and  (17),  one  has  for  m=2: 


IX/s  |X2|R  |X3|rt  |X4 

ixirtix2p,!x3rsix4 


7  |x5|rs  |X6|r«,  |x7p,  |X8|r*,  IX/t  |x10 
id  ix5rt  ix6p/  ix7r/  ix8p>,  |x9d  ixI0| 


(39) 

(40) 


The  last  three  columns  of  Table  1  list  the  errors  in  the  estimates  for  BP  and 
BV for  the  three  sets  of  times,  [TO],  [Tl],  [T2],  In  the  case  m=l,  the  esti¬ 
mates  are  given  for  the  x-component  of  BP  for  B1,Bl1,B2g.  The  maximum 
error  for  the  latter,  over  [T2],  is  2533  meters,  occurring  at  an  IM  intercept 
time  of  TE(82)=1715.  The  results  are  considered  satisfactory  considering  the 
complex  nature  of  the  functions  (see  Figure  1A  below).  In  addition,  using 
the  interpolated  values,  at  TE(82)=1715,  for  BP  and  BV,  the  miss  distance 
of  the  IM  at  the  target  is  100.24  meters,  which  should  also  be  acceptable. 

The  figures  that  follow,  1A,  IB,  1C,  ...,  6A,  6B,  6C,  show  the  six  different 
types  of  functions  that  were  fitted  and  the  errors  that  were  observed.  Figures 
mA,  m=l,2,3,4,5,6,  contain  plots  of  the  positive  term  in  (36)  as  a  function 
of  t  e  [T2],  Figures  mB  show  the  plots  of  the  sum  of  the  two  negative  terms 
in  (36)  taken  with  plus  signs.  Figures  mC  contain  the  errors  observed,  (36). 
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,6  m  =  1 ,  Jrm=  (2  9  4.85  0.9  2.275) 


Figure  1A.  x-Component  of  Input  Burnout  Position  versus  Time  [T2] 


x1Q6  m  =  1 ,  jrm=  (2  9  4.85  0.9  2.275) 


Figure  IB.  x-Component  of  Output  Burnout  Position  versus  Time  [T2] 


li 


ErrorT2 
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m  =  1 ,  Jrm=  (2  9  4.85  0.9  2.275) 
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B65=  2,  Exp  values  =  5  7  2.275  2.45  1 


Figure  2A.  y-Component  of  Input  Burnout  Position  versus  Time  [T2] 


B65=  2,  Exp  values  =  5  7  2.275  2.45  1 


Figure  2B.  y-Component  of  Output  Burnout  Position  versus  Time  [T2] 
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m  =  2,  Jrm=  (5  7  2.275  2.45  1) 


T  1  1  1 

*  Oa.ICUIa.IGQ  tiTOl  1  ^ 

1000  -  1 


-1000 


900  1100  1300  1500  1700  1900  2100  2300  2500 

Time  (seconds) 

Figure  2C.  Error  in  y-Component  of  Burnout  Position  versus  Time  [T2] 
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Burnout  position,  z-component  (meters)  ^  Burnout  position,  z-component  (meters) 
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,6  m  =  3,  Jrm=  (4  8  1.625  0.9  1) 


igure  3A.  z-Component  of  Input  Burnout  Position  versus  Time  [T2] 


x1Q6  m  =  3,  Jrm=  (4  8  1.625  0.9  1) 
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m  =  3,  Jrm  =  (4  8  1.625  0.9  1) 


Figure  3C.  Error  in  z-Component  of  Burnout  Position  versus  Time  [T2] 
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m  =  4,  Jrm  =  (2  7  3  1.5  0.725) 


Figure  4A.  x-Component  of  Input  Burnout  Velocity  versus  Time  [T2] 


m  =  4,  jrm  =  (2  7  3  1.5  0.725) 
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m  =  4,  Jrm  =  (2 


3  1.5  0.725) 


Figure  4C.  Error  in  x-Component  of  Burnout  Velocity  versus  Time  [T2] 
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m  =  5,  Jrm=  (2  6  3.825  1.4  1) 


Figure  5A.  y-Component  of  Input  Burnout  Velocity  versus  Time  [T2] 


m  =  5,  jrm=  (2  6  3.825  1.4  1) 
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m  =  5,  Jrm  =  (2  6  3.825  1.4  1) 


1  1  \ 

- * - Calculated  ErrorT2 

- 

0.4 


900  1100  1300  1500  1700  1900  2100  2300  2500 

Time  (seconds) 

Figure  5C.  Error  in  y-Component  of  Burnout  Velocity  versus  Time  [T2] 
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m  =  6,  Jrm=  (4  5  2.7  1.625  1) 


igure  6A.  z-Component  of  Input  Burnout  Velocity  versus  Time  [T2] 


m  =  6,  jrm=  (4  5  2.7  1.625  1) 
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m  =  6,  Jrm=  (4  5  2.7  1.625  1) 


Figure  6C.  Error  in  z-Component  of  Burnout  Velocity  versus  Time  [T2] 
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ATHENS  GA  30602 

DR  WILLIAM  GRIFFITH  1 

DEPT  OF  STATISTICS 
UNIVERSITY  OF  KENTUCKY 
LEXINGTON  KY  40506-0027 

DR  REBECCA  W  DOERGE  1 

DEPT  OF  STATISTICS 

PURDUE  UNIVERSITY 

1399  MATHEMATICAL  SCIENCE  BLDG 

WEST  LAFAYETTE  IN  47907 

DR  LARRY  WASSERMAN  1 

DEPT  OF  STATISTICS 
CARNEGIE  MELLON  UNIVERSITY 
PITTSBURGH  PA  15213 

DR  SATISH  IYENGAR  1 

DEPT  OF  STATISTICS 
UNIVERSITY  OF  PITTSBURGH 
2703  CATHEDRAL  OF  LEARNING 
PITTSBURGH  PA  15260 

DR  DIPAK  K  DEY  1 

DEPT  OF  STATISTICS 

UNIVERSITY  OF  CONNECTICUT 

196  AUDITORIUM  ROAD 

STORRS  CT  06269 

DR  ANDREW  BARRON  1 

DEPT  OF  STATISTICS 

YALE  UNIVERSITY 

PO  BOX  208290 

NEW  HAVEN  CT  06520-8290 


Copies 

DR  KESAR  SINGH  1 

DEPT  OF  STATISTICS 

RUTGERS  UNIVERSITY  OF  NEW  JERSEY 
501  HILL  CENTER  BUSCH  CAMPUS 
PISCATAWAY  NJ  08854 

DR  ROBERT  SMYTHE  1 

DEPT  OF  STATISTICS 
GEORGE  WASHINGTON  UNIVERSITY 
WASHINGTON  DC  20052 

DR  STEPHEN  L  BIEBER  1 

DEPT  OF  STATISTICS 

UNIVERSITY  OF  WYOMING 

PO  BOX  3332 

LARAMIE  WY  82071-3332 

DR  HOWARD  B  CHRISTENSEN  1 

DEPT  OF  STATISTICS 

BRIGHAM  YOUNG  UNIVERSITY 

230A  TMCB 

PROVO  UT  84602 

DR  ELGIN  PERRY  1 

2000  KING’S  LANDING  RD 
HUNTINGTOWN  MD  20639-9743 

JOHN  T  CARR  1 

VIRGINIA  TECH  CORPORATE  RESEARCH  CENTER 


1800  KRAFT  DRIVE  SUITE  111 
BLACKSBURG  VA  24060 

INTERNAL 

CX7  TECHNICAL  LIBRARY  3 

G30  CHARLES  GARNETT  1 

G72  NGA  PHAM  1 

Q31  DR  JOHN  CRIGLER  1 

Q32  DR  MICHAEL  RUDZINSKY  1 

Q33  DR  JEFFREY  SOLKA  1 

W10  EMMANUEL  SKAMANGAS  1 

Wll  STEVEN  ANDERSON  1 

Wll  WILLIAM  ORMSBY  1 

W14  DAVID  CLAWSON  1 

W20  DAVID  RICHARDSON  1 

W25  JOHN  BECK  1 

W25  DR  ARMIDO  DIDONATO  2 

W25  DR  ROBERT  MCDEVITT  1 

W25  DR  DEMETRIOS  SERAKOS  1 

W25  DANA  HUMPHREY  1 

W25  DR  MICHAEL  HABGOOD  1 

W60  ROBERT  HILL  1 

W60  SAID  SAADI  1 
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